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ABSTRACT 

We have discovered five millisecond pulsars (MSPs) in a survey of 14 unidentified Fermi-LAT sources 
in the southern sky using the Parkes radio telescope. PSRs J0101-6422, J1514-4946, and J1902-5105 
reside in binaries, while PSRs J1658-5324 and J1747-4036 are isolated. Using an ephemeris derived 
from timing observations of PSR J0101-6422 (P=2.57ms, DM=12pccm -3 ), we have detected 7-ray 
pulsations and measured its proper motion. Its 7-ray spectrum (a power law of T = 0.9 with a cutoff 
at 1.6 GeV) and efficiency are typical of other MSPs, but its radio and 7-ray light curves challenge 
simple geometric models of emission. The high success rate of this survey — enabled by selecting 7- 
ray sources based on their detailed spectral characteristics — and other similarly successful searches 
indicate that a substantial fraction of the local population of MSPs may soon be known. 

Subject headings: gamma rays: general — pulsars: individual (PSR J0101-6422) 



1. INTRODUCTION 



The Large Area Telescope (LAT, lAtwood et aHl2009f ) 
of the Fermi Gamma-ray Space Telescope is a pulsar 
detect or par excellence. Supported by a radio cam- 
paign (|Smith et all [20081 providing ephemerides for co- 
herent folding of LAT photons, the LAT has identified 
7-ra y pulsations from m any normal (unrecycled) pulsars 
(e.g. lAbdo et al.ll2010cf ) and millisecond pulsars (MSPs, 
lAbdo et all l2009af ). Additionally, tens of new pulsars 
have been discovered i n "blind" searches for periodicity 
in th e LAT data (e.g. lAbdo et al.ll2009bl: iPletsch et at] 
[201211 . 

A third, indirect method of detection has been ex- 
tremely successful in discovering new pulsars, espe- 
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cially MSPO Gamma -ray pulsars — both young and 
recycled — are stable emitters and nearly all have a char- 
acteristic spectrum dN/dE oc E~ T exp(— E/E c ) wit h 
T < 2 and 1 < E c /GeV < 10 (|Abdo et al.l I2010cl) . 
making unidentified, nonvariable LAT sources with such 
spectra good pulsar candidates. The fine angular resolu- 
tion and large effective area of the LAT allow a typical 
source localization of ~10', about the 1 GHz beam size of 
100 m class radio telescopes. This happy coincidence en- 
ables deep pulsation searches with a single pointing. Dis- 
covery of a pulsar and a subsequent timing campaign can 
"close the loop" by providing an ephemeris with which to 
fold LAT photons and resolve 7-ray pulsations. Indeed, 
since most MSPs reside in binaries, initial detection and 
characterization of the orbit at longer wavelengths is the 
only way 7-ray pulsations can be detected from these 
pulsars. 

Due to sensitivity limitations, LAT-guided surveys nat- 
urally target the nearly-isotropic population of nearby 
MSPs. The deep radio exposures afforded by the effi- 
cient target selection offer the tantalizing possibility of 
completely cataloguing the local field MSPs whose radio 
and 7-ray beams cross the Earth. This relatively com- 
plete sample will be valuable in constraining the emission 
mechanisms of MSPs and in characterizing the evolution 
of their binary progenitors. 

Towards this end, we performed a small, targeted sur- 
vey at the CSIRO Parkes t elescope, compleme ntary to a 
Parkes search described in IKeith et al.l (|201 lh that had 
different source selection criteria. 

2. TARGET SELECTION 

We began with the set of southern (Deck < —40°; 
northern sources are visible to the more sensitive Green 
Bank Telescope), nonvariable sources in a preliminary 
version of the 1FGL catalog of 7-ray sources (jAbdo et al.l 

17 See https:/ /confluence. slac. Stanford. edu/display/GLAMCOG/Public+List 
Detected+Gamma-Ray+Pulsars for an up-to-date list of LAT- 
detected pulsars 
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I2010H) . We further restricted consideration to sources 
whose LAT position estimates had a 95% confidence 
error radius <7', the HWHM of the 1.4 GHz beam of 
the Parkes telescope, and which had no plausible blazar 
counterpart. We classified the remaining sources by vi- 
sual inspection of the 7-ray spectrum. Sources with spec- 
tral shapes resembling known pulsars (see above) were se- 
lected for observation, while those sources best described 
by power laws or with significant spectral breaks but 
T > 2 were deemed likely to be blazars and were not 
selected. Due to the difficulty of spectral modeling and 
localization of LAT sources in the Galactic plane, we dis- 
carded sources with |6| < 5°. Our final list comprised 14 
good candidates. The source positions (see Table 1) were 
computed from the same analysis used to estimate the 
spectral shape, but in all cases the difference in position 
with 1FGL is well within both the 1FGL 95% error con- 
tour and the Parkes beam. 

3. RADIO SEARCHES 

We observed the 14 selected sources (S}2]) at the Parkes 
64-m radio telescope between 2009 November 25 and 
December 8. Each source was observed for between 
1 and 2 hours (see Table 1) using the center beam 
of the multibeam receiver. We recorded 1-bit-digitized 
total-power samples every 125 /is from each of 512 fre- 
quency channels spanning 256 MHz of band centered on 
139 MHz, writing the data to disk for off-line analysis 
(see iManchest er et all (|2001[ ) for more details on the re- 
ceiver and data-acquisition system.) 

The data were analyzed on the koala computer 
cluster at Columbia University using standard pulsar 
search techniqu es implemented in the PRESTO package 
(Ransom 2001). We dedispersed each data set ideally up 
to a dispersion measure of DM = 270pccm~ 3 , and when 
searching for periodic signals we allowed for accelerated 
signals caused by pulsar motion in a binary system. This 
was parameterized within PRESTO by the parameter 
zmax = 50, which is related to the maximum number 
of bins a signal can drift within the Fourier spectrum 
and still be properly detected by the search method. All 
of the data were analyzed in this manner within about a 
week of collection, and on 2009 December 4 we found our 
first MSP. Within a few more days we confirmed it and 
four other new MSPs. A 6th MSP was detected unbias- 
ed ly in our sample, bu t it had been previously discovered 
bv lKeith et all (pill see Table 1). 

In a few cases, the pulsars were not confirmed on the 
first attempt and required subsequent observations, ow- 
ing to their faintness and scintillation in the interstellar 
medium. We reanalyzed the data by increasing the ac- 
celeration search space to zmax = 200 and by analyzing 
only the first 30 minutes of each data file. The latter 
approach improves our sensitivity to highly accelerated 
pulsars (e.g., for an MSP in a 1-day binary, a 2-hour ob- 
servation cannot be corrected ideally under the assump- 
tion of a constant acceleration, no matter how large the 
zmax used), but no new pulsars were detected. 

The nominal sensitivity of our searches, for a 2 hr ob- 
servation and assuming a pulsar duty cycle of 20%, is 
0.11 mJy for P >2ms and DM < 40pccm~ 3 , degrading 
gradually for shorter periods and larger DMs. While in 
detail this depends on the sky background temperature, 
the contribution from the Galaxy at our observing fre- 



quency only differs by a maximum of 3 K for different 
locations, out of a total system equivalent temperature 
of nearly 30 K, so that we provide an average flux density 
limit. Two of the non-detections have limits of 0.12 mJy 
and 0.14 mJy owing to shorter integration times (see Ta- 
ble [lj . Because of scintillation, a pulsar with a larger 
average flux than our nominal limits might not be de- 
tectable in a particular observation, and vice versa. In 
any case, the nominal flux density represents the most 
sensitive searches for MSPs done at Parkes, especially 
considering that these LAT-selected targets are largely 
expected to be relatively nearby. 

Characterization of four of the MSPs requires a longer 
radio timing campaign and/or LAT dataset than avail- 
able for this work. We defer discussion of these MSPs 
(and of scintillation effects on sensitivity, which like- 
wise requires repeated observations) to Camilo et al. (in 
preparation). The remaining pulsar presents a peculiar 
light curve and we describe it further below. 

4. PSR J0101-6422 
4.1. Timing Solution 

After the discovery observations on 2009 November 25 
and confirmation on December 8, we began regular tim- 
ing observations of PSR J0101-6422 at Parkes. We ob- 
served the pulsar on 35 days through 2011 November 9, 
detecting it on 28 occasions; on the remaining 7 days it 
was too faint to detect, due to interstellar scintillation 
(DM=12 pc cm -3 .) Each observation was typically one 
hour, with the same receiver and data acquisition system 
used in the search observations. Using TEMP023 with 
the 28 times-of- arrival, we obtained a phase-connected 
timing solution whose pa rameters are giv en in Table [5J 

The Shklovskii effect (|Shklovskiil 11970ft increases the 
measured value of P = 5.2 x 10~ 21 beyond that intrin- 
sic to the pulsar by AP = Pv\/dc, with v± the source 
velocity transverse to the line-of-sight. For nearby ob- 
jects, it may dominate the true spindown ra te and lead to 
overe stimates of the spindown luminosity ([Camilo et al.l 
1994). The observed proper motion, at the 0.55 kpc DM 
distance ()Cordes fc Laziol |2()()2: , implies a true Pi — 
(4.3 ± 0.1) x 10~ 21 . At this distance, the indicated 
v 1 = 41 kms -1 is typical for an MSP (see lNice" fc Taylor! 
fl995h . 

4.2. Radio Profile and Polarimetry 

We made polarimetric and flux-calibrated observations 
of PSR J0101-6422, using the PDFB3 digit al filterbank 
The data were analyzed with PSRCHIVE (jHotan et al.l 
120041 ). The best full-Stokes profile that we obtained is 
shown in Figure [T] (the flux density of this profile is not 
necessarily representative of the average pulsar intensity, 
which varies greatly owing to interstellar scintillation). 
The main pulse is linearly polarized at the ^15% level, 
but the low signal-to-noise ratio prevents a useful de- 
termination of rotation measure (the best-fit value is 
~ +10radm~ 2 , but within the uncertainties is consis- 
tent with zero). These data show that the subsidiary 
pulse component is composed of two outer peaks joined 
by a low bridge of radio emission. 

18 http:/ /www. atnf.csiro.au/research/pulsar/tempo2 
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Table 1 

Results of Radio Searches of 14 1FGL Sources at Parkes 



11 GL Name 


R.A. a 


Decl. a 


1 


b 


r 


E c 


Obs. Time 


Period 


DM 


Distance 


Binary 




(J20U0.U) 


(J2000.0J 


(deg) 


(deg) 




(GeV) 


(hr) 


(ms) 


(pc cm -3 ) 


(kpc) 




J0101. 0-6423 


01 h 00 m 58?l 


-64°24'03" 


301.2 


-52.7 


1.3 


2.3 


1.0 


2.57 


11.93 


0.55 


yes 


J0603. 0-4012 


06 h 03 m 04=9 


-40°11'02" 


246.8 


-25.9 


1.9 


6.5 


2.0 










J0933.9-6228 


09 h 33 m 58=5 


-62° 27' 54" 


282.2 


-7.8 


0.1 


1.2 


2.0 














J1036.2-6719 


10 h 36 m 16?9 


-67°20'34" 


290.4 


-7.8 


1.5 


2.2 


2.0 


— 


— 


— 


— 


J1227.9-4852 


12 h 27 m 50?5 


-48°51'54" 


298.9 


13.8 


1.8 


1.6 


2.0 










J1232.2-5118 


X2 h 31 m 49H 


-51° 18' 50" 


299.8 


11.4 


1.8 


2.3 


2.0 










J1514.1-4945 


15 h 14 m 05?7 


-49°45'32" 


325.2 


6.8 


1.7 


5.9 


2.0 


3.59 


31.5 


0.9 


yes 


J1624.0-4041 


16 h 24 m 06?2 


-40°40'48" 


340.6 


6.2 


2.1 


3.3 


2.0 










J1658.8-5317 


16 h 58 m 43?2 


-53° 17' 45" 


335.0 


-6.6 


2.1 


1.8 


1.3 


2.44 


30.9 


0.9 


no 


J1743.8-7620 


17 h 43 m 44?6 


-76°20'42" 


317.1 


-22.5 


1.2 


2.0 


1.5 










J1747.4-4035 


17 h 47 m 29 B 1 


-40°36'07" 


350.2 


-6.4 


1.5 


3.4 


1.4 


1.65 


152.9 


3.3 


no 


J1902.0-5110 


19 h 02 m 05?5 


-51°09'43" 


345.6 


-22.4 


1.7 


4.4 


1.2 


1.74 


36.3 


1.2 


yes 


J2039.4-5621 


20 h 39 m 30?5 


-56°20'42" 


341.2 


-37.1 


1.6 


2.7 


1.2 










J2241.9-5236 c 


22 h 41 m 52?4 


-52°37'37" 


337.4 


-54.9 


1.6 


3.6 


2.0 


2.19 


11 


0.5 


yes 



a Parkes telescope pointing position. 

b Dist ances computed f rom measured DM using the NE2001 model (Cordcs & Lazio 2002). 
c See lKeith et al l lIMTH l. 
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Figure 1. Polarimetric pulse profile of PSR J0101-6422 at 
1.4 GHz, displayed with 128 phase bins, based on 6hr of Parkes 
PDFB3 data using 256 MHz of bandwidth. The black trace corre- 
sponds to total intensity, red to linear polarization, blue to circular. 
In the upper window, the position angle of linear polarization is 
plotted for bins with linear signal- to- noise ratio > 3. The pulse 
peak is displayed with arbitrary phase, and the mean flux density 
is ~ 0.2 mjy. 

4.3. Gamma-ray Profile 

To characterize the 7-ray profile, we selected photons 
collected between 2008 Aug 4 and 2011 Aug 1 with 
reconstructed energies 0.3-3 GeV lying within 1.1° of 
the timing position. This selection balances complete- 
ness (including as many pulsed photons as possible) and 
pulsed signal-to-nois e ratio as measured by the H-tcst 
(|de Jager et al.|[l989l) . We apply the same data process- 
ing as described in M4.41 omitting the horizon cut to in- 
crease the livetime by about 20%, at the expense of a 
slightly increased background. 

The light curve corresponding to this extraction ap- 
pears in Figure [2] If 7 rays are produced above the null 
charge surface (the locus of points where the magnetic 
field is orthogonal to the pulsar spin axis), they appear 
to an observer in the opposite hemisphere, whereas ra- 
dio emission from low altitudes is beamed into the same 
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Figure 2. Phase-aligned 7-ray (blue histogram) and radio (red 
trace) pulse profiles of PSR J0101-6422. The radio profile is 
summed from 28 timing detections, while the 7-ray profile cor- 
responds to the optimal aperture described in the main text. The 
background estimation (dashed horizontal li ne) is derived by count- 
ing photons in the off-pulse phase window ( c!4.4ll and is consistent 
with the off-pulse background level obtained via spectral analysis. 
The vertical dashed line indicates the phase (</> = 0.5) from which 
the offset of "PI" is measured. 

hemisphere (see iRomani fc Yadigaroglul 119951 ). In this 
scenario, the line-of-sight to PSR J0101-6422 would in- 
tercept a cone of radio emission in one hemisphere to 
produce the weaker radio peak at <fi ~ 0.5 and a cone of 
radio and 7 emission in the other to produce the two 7 
peaks and bright radio peak at <f> ~ 1. In keeping with 
this picture and the convention established in the litera- 
ture, we identify the 7 peak at <f> ~ 0.8 as "PI" and that 
at <£~0.1 as "P2". 

The relative phas ing of the peaks is of interest (e.g. 
iWatters et all 120091 ) . The radio-to-7 offset, 5, is mea- 
sured from the center of the leading radio peak (taken as 
(p = 0.5, see Figure [5]). We estimated the phase of the 
7-ray peaks using unbinned maximum likelihood with 
a two-sided Lorentzian model for the peaks. These esti- 
mates yield 5 = 0.33±0.01 s t a t±0.03 sys t and a 7-ray peak- 
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kpc. Best-fit values of /q from the TPC and OG models 



5.ip are 0.7 and 0.9. This 7-ra y efficiency is in ac cord 



with other LAT-detected MSPs fAbdo et al.ll2009al) 

4.5. X-ray and Optical Observations 

A 3.3 ks Swift observation with the X-ray Telescope in 
PC mode was obtained on 2009 Nov 21 to search for an 
X-ray counterpart to the then-unidentified 1FGL source. 
No significant counterpart is detected, with a 3ct upper 
limit at the pulsar position of 2.0 x 10 _3 ctss _1 (0.5-8 
keV). Assuming a power-law spectrum with photon index 
T = 1.5 for a column with A H = 3.7 x 10 20 cm- 2 (10 
HI atoms per free electron along the line-of-sight) yields 
an unabsorbed flux limit of 1.1 x 10~ 13 erg cm~ 2 s -1 , or 
efficiency Lx/E < 5.0 x 10~ 4 at 0.55 kpc. This limit is 
com parable to the detec ted X-ray flux from other MSPs 
(e.g. lMarelli et al.l[20Tl . 

We searched four DSS plates and found no optical 
counterpart consistent with the pulsar's position, sug- 
gesting to > 21. 

5. DISCUSSION 

5.1. PSR J0101-6422 

At first blush, PSR J0101-6422 is an unremarkable 
member of the growing population of 7-ray MSPs. The 
orbital characteristics suggest an epoch of Roche lobe 
overflow from an evolved 1-2 Mq companion which later 
became the ~0.2 He white dwarf curr ently inferred 
in the system (e.g. lTauris &; Savonijdligggt ). However, as 
we discuss below, PSR J0101-6422's light curve is chal- 
lenging to explain with simple geometric models that tie 
the radio and 7 emission to particular regions of the mag- 
netosphere. 

Remarkably, a simple picture of 7-ray emission aris- 
ing in the outer magnetosphere ap pears to describe the 
majority of both young (th ough seelRomani 
and recycled pulsars (e.g. iVenter et al. 2009). On the 
other hand, while radio emission in young pulsars and 
some MSPs appears to come from lower altitudes, sev- 
eral MSPs have p hase-aligned 7 and radio peaks (e.g. 
lAbdo et~al . 2010a), implying a joint emission site at high 
altitude. These observations raise the prospect of MSPs 
whose radio emission combines "traditional" polar cap 
emission with a higher-altitude component. 

To determine which scenario best describes the emis- 
sion pattern of PSR J0101-6422, we performed joint 
fit s to the r adio a nd 7-ray data following the method 
of Uohnsonl (|2011f ) . We modelled the rad io as a sim- 
ple "cone" component (|Storv et al.l 120071 ) and the 7 
ra ys using both the "tw o-pole caustic" (TPC) model 
of lDvks & Rudakl (|2003l N.B. we adopt a lar ger max- 
imum cylindrical radius of 0.95 of the light cylinder 
and the "outer gap" (OG) model (jCheng et all 11986 
iRomani fc Yadigaroglul TT995T ) . The angle between the 
pulsar spin axis and the line-of-sight, £, is unknown a 
priori, but if the radio emission arises from low altitudes, 
then the presence of two peaks separated by roughly half 
of a rotation indicates the inclination of the magnetic 
axis from the spin axis (a) cannot be small. The best-fit 
light curves for the two 7-ray models appear in Figure [31 
The low a = 26° preferred for the TPC model cannot 
produce the radio interpulse and produces too much off- 
pWse 7-ray emission. The OG model, with a best fit at 



0.8 1 1.2 
Rotation Phase 

Figure 3. The best-fit model light curves for the two-pole caus- 
tic (TPC) and outer gap (OG) models (see fl57Tl ) The best fit 
geometries («,£) for TPC and OG are (26°, 79°) and (90°, 36°). 

to-peak separation of A = 0.26 ± 0.03 sta t ± 0.02 syst . The 
DM-induced uncertainty in S is < 0.001P. The system- 
atic uncertainty for both 5 and A is estimated from the 
scatter obtained by fitting alternative functional forms 
for the peaks. 

4.4. Gamma-ray Spectrum 

To measure the spectrum of psr, we analyzed low back- 
ground DIFFUSE class "Pass 6" (jAtwood et al.l2009tl LAT 
data collected between 2008 Aug 4 and 2011 Apr 17. We 
filtered periods where the observatory's rocking angle ex- 
ceeded 52° or the Earth's limb impinged upon the field 
of view (horizon cut, requires zenith angle < 100°). We 
modelled the sensitivity to these events with the flight- 
corrected P6_V11_DIFFUSE instrument response function. 
To compute background contributions, we used a prelim- 
inary version of the 2FGL catalog and its accompany- 
ing diffuse emission modelJ^I. We selected events within 
10° and re-fit the spectra of point sources within 8° and 
the norm alizations o f the diffuse models. Fits employed 
pointlike (jKerrl 120111 ). a binned maximum likelihood al- 
gorithm. 

The phase-averaged spectrum of PSR J0101-6422 is 
well described by an exponentially cutoff power law, 
dN/dE = N {E/E )~ r exp(-E/E c ), whose parame- 
ters appear in Table [2J We also considered separately 
the spectra of the three major components, viz. P2 
(0 < (j) < 0.33), the off-pulse (0.33 < 4> < 0.63), and 
PI (0.63 < 4> < 1). The parameter uncertainties are 
relatively large, and accordingly we found no significant 
difference in the spectral shapes of the two peaks. Fur- 
ther, we detected no significant emission in the off-pulse 
phase window and derived a 95% confidence upper limit 
on the energy flux of a point source with a power law 
spectrum with T = 2.2, obtaining 6 x 10 -12 erg cm~ 2 
s -1 (> 100 MeV; scaled up by 3.3 to the full phase win- 
dow.) 

The observed flux is related to the 7-ray luminosity 
by F 7 = L 7 /47r /n d 2 , where /q is a n unknown beaming 
factor (typically of order unity, e.g. iWatters et al.|[2009l) 
implying a phase-averaged, isotropic (/n = 1) 7-ray lu- 
minosity L 1 = 4.0 x 10 32 ergs _:L = 0.04 E at d = 0.55 

19 http: / /fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels 



5 Millisecond Pulsars in Fermi Sources 



Table 2 

Measured and Derived Parameters for PSR J0101-6422 

Parameter Value 

Right ascension, R.A. (J2000.0) 01 h 0l m lif 1163(2) 

Declination, decl. (J2000.0) -64°22'30'.'171(2) 

Proper motion in R.A. X cos decl. (mas yr ) .... 10(1) 

Proper motion in decl. (mas yr _1 ) —12(2) 

Galactic longitude, I (deg.) 301.19 

Galactic latitude, b (deg.) -52.72 

Spin period, P (ms) 2.5731519721683(2) 

Period derivative (apparent), P 5.16(3) X 10~ 21 

Epoch of period (MJD) 55520 

Orbital period, P b (days) 1.787596706(2) 

Epoch of periastron (MJD) 55162.4011764(3) 

Semi-major axis (It s) 1.701046(2) 

Eccentricity < 1.5 X 10~ 5 

Dispersion measure, DM (peem" 3 ) 11.926(1) 

Spin-down luminosity a , E (ergs -1 ) 1.0 X 10 34 -16% 

Characteristic age a , t c (yr) 9 X 10 9 +20% 

Surface dipole magnetic field strength (Gauss) a . . . 1.1 X 10 7 —9% 

Mean flux density 13 at 1.4 GHz, Si. 4 (mjy) 0.28 ± 0.06 

Radio-7-ray profile offset, S (P) 0.35 ± 0.01 ± 0.05 

7-ray profile peak-to-peak separation, A (P) 0.26 ± 0.03 ± 0.02 

7-ray (> 0.2 GeV) photon index, T 0.9 ± 0.3 ± 0.2 

7-ray cut-off energy, E c (GeV) 1.6 ± 0.4 ± 0.3 

Photon flux (> 0.1 GeV) (10 -9 cm" 2 s" 1 ) 9.5 ± l.sj^' 

Energy flux (> 0.1 GeV), F 7 (10 -11 erg cm" 2 s" 1 ) 1.1 ±0.1 ± 0.1 

Note. — Timing solution parameters given in TDB relative to the DE405 plane- 
tary ephemeris. Numbers in parentheses represent the measured 1 a TEMP02 timing 
uncertainties on the last digits quoted. For 7-ray parameters, the first uncertainty is 
statistical and the second systematic. 

a Corrected for Shklovskii effect at 0.55 kpc; the percentage change from the nominal 
value is given. 

b Average flux for 28 timing detections. 

a = 90° (orthogonal rotator), produces two radio peaks | b| > 5°, their efficiency jumps to 3/8. iKeith et"ai"1 

but fails to produce two clear 7 peaks and the radio in- (|2011[ ) employed some spectral information and detected 

terpulse morphology. MSPs unbiasedly in 1/1 f pointings (1/4 for \b\ > 5°). 

Since neither model faithfully reproduces the observed A short list of 7 spectrally-exceptional candidates drawn 

light curves, we propose that one or both of the radio up for a Nangay radio telescope survey included 5 MSPs, 

peaks may originate at higher altitude. Emission over a three of which were discovered from Nancay o bserva- 

range of altitudes would also explain the lo w observed tions flCognard et alJl20lH iGuillemot et aT1l2011[ ). And, 

linea r polarization of PSR JOf 01-6422 (e.g. iDvks et"aTI iHessels et al.l (|2011l ) used nearly identical selection crite- 

2004) and admit smaller a values. Alternatively, the fail- ria as this work in a Green Bank Telescope survey at 350 

ure of geometric models may indicate that the physical MHz and found MSPs in 13/49 pointings, all but one at 

details of the MSP magnetosphere — e.g., currents, mul- |6| > 5°. the 

tipolar fields, and plasma loading of the field lines — and Targeting pulsar- like sources is clearly efficient. How- 

the radio and 7 emission mechanisms play an important ever, new LAT sources will lie at the sensitivity threshold 

role. and the limited statistics will not admit classification of 

spectral shape and variability. Nonetheless, we argue 

5.2. Survey Implications that future radio searches should omit detailed ranking 

„,.,, nM» ■ a- ,,. • ,. , ,,n n i . schemes and target any LAT source off the Galactic plane 

With 6/14 pointings resulting m unbiased MSP detcc- , , , , . , , .,, , v 

.' ° r ., r , 1 . 1 • 1 • 1 rr> • that has not been associated with a known blazar. 

tions, this survey was fruitful, though tins high efficiency TTr , ,. , , . , , ,. . , ,. , , , 

u i 1 T j j • i.- • • i j We beheve this substantiafiy increased target list (and 

may be luck. Indeed, m our timing campaign, we only de- , , ,. \ ....•«•> 1 T .1 ± i 

1 inon Tmm no roc i-L x j 1 telesco pe time) is lustified. In the 21 GL source catalog 

tected PSR J0101-6422 on 28 of 35 attempts, and severaf \ A i , r , — ni onn h ^ ocn • * j 

, ,. • 1 • 1 1 1 cj.1? 1 • (lAbdo et al. 20111), there are ^350 unassociatcd sources 

observations required prior knowledge of the ephemeris ' ~r PS — a A. n 1 j.- 1 ..cti , c 

r , , , . ^ , ^ c . j.- j. more than 5 on the Gaiactic pfane. MSPs account for 

tor detection. Ihus, because of scmtiffation, we estimate -n? <• j T »rp 1 <■ on nor, 

that PSR J0101-6422 is only detectable in 2/3 observa- 5% °f as K S0C f ted , ^, SOurces ' 80 15 7 2 ? ME- 8 

n ,. , , . , ,.' , , may remain to be found m these unassociatcd 21 GL 

tions, depending somewhat on integration time. Many J n,.^ , „_„ ... .. ,.. , 

c u irc-n ■ l- u- • r. it. i • . • v • sources. With only 350 positions to monitor, multiple 

held MSPs arc m eclipsing binaries, further diminishing , , J ,. r „ , 

, r j . . • • • 1 i_ deep observations can ameliorate the confounding effects 

the probability of detection in a smgfe observation. c ■ , . ln , . , ,. , ., . . c °. . c 

TX f , r . ■ i5 of scmtiffation and echpses whife requiring a traction of 

We befieve the strongest factor m our search s success , , , . ? f, 1j.n1 o 

, ,.° ., . y-, rr: — n the time of a comparabfy compfete aff-sky survey. By 

was the source seiection criteria. E.g., I Ransom et af.l . ,f . , . . n 1 i^-u 

.1 , 1,1 • 1 1 x Arn — n — T~ '■ — 1 — 'r comparison, the otherwise prodigious Parkes muftibeam 

2011) seiected nonvariabfe LAI sources (but mcfud- r r , « , „ 1f7 . An nnn . ,. 

' , 1 . c \ j C j n/ron • survey found onfy 17 MSPs m over 40,000 pointings 

mg no spectraf shape information) and found MSPs in , tk — rr^ , ' i ii onn/ih t jj-j-- l i.v. ■ • • 

o /ot • j.- r> : i_ . i , , , (e.g. laufkner et af. 2004 . In addition to the intrinsic 

3/25 pointings. But when considering onfy targets with v ' 
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interest in discovering new MSPs, by thoroughly search- 
ing all such LAT sources, we will gain confidence that the 
MSPs we have found can truly provide a volume-limited 
sample of the energetic, 7-loud MSPs, an invaluable step 
in the study of the MSP population of the Milky Way. 
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